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Abstract. V819 Herculis is a well-studied triple star system consisting of a "wide" pair with 5.5 year period, one component 
of which is a 2.2-day period eclipsing single-line spectroscopic binary. Differential astrometry measurements from the Palomar 
High-precision Astrometric Search for Exoplanet Systems (PHASES) are presented and used to determine a relative inclination 
between the short- and long-period orbits of 23.6 ± 4.9 degrees. This represents only the sixth unambiguous determination of 
the mutual inclination of orbits in a hierarchical triple system. This result is combined with those for the other five systems 
for analysis of the observed distribution of mutual inclinations in nearby triple systems. It is found that this distribution is 
different than that which one would expect from random orientations with statistical significance at the 94% level; implications 
for studying the spatial distribution of angular momentum in star forming regions is discussed. 
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1. Introduction 

Determinations of the mutual inclinations of the two orbits in 
hierarchical triple stellar systems are rare, with previously only 
four unambiguous determinations available in the liter ature 
iLestrade et al.Lll993llHein^ll996llHurnmel et allEool . but 
valuable; the dynamical relaxation process undergone by mul- 
tiples after formation is expected to le ave a statistical "finger- 
print" in the distribution of inclinations | Ste rzik and Tokovinini 
2002]. In addition, well-characterized stellar multiples repre- 
sent excellent opportunities to test and challenge stellar models 
under stringent constraints of common age and metallicity. 

V819 Herculis (HR 6469, HD 157482; V = 5.6, K = 4.1) 
is a triple system consisting of an evolved star (G7 III-IV; this 
appears brighter in V and will be referred to as the A com- 
ponent) in an eccentric 5.5 year orbit together with a close 
(P„ * 2.2 days) pair of main sequence F stars. The close 
pair (Ba and Bb) is in an edge-on orbit and exhibits shallow 
eclipses. A combination of radial velocity, speckle interferom- 
etry, eclipse timing and light-curve fitting has made it possible 
to accurately determine most of the interesting system param- 
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eters, including masses, radii and distance IScarfe et ail 1 1994 
Ivan Hamme et all 1 19941 IWasson et all Il994ll with accuracies 
of a few percent. However, until now it has not been possi- 
ble to determine the mutual inclination of the orbits. The V819 
Her system is listed as a chrom osphericallv active binary in 
the catalog by IStrassmei er et al ] il993ll : it exhibits Ca H and 
K em ission and has b een detected in X-r ays IPempsev et all 
Il993ll but not in radio |Drake et al., 1989]. The system is vari- 
able with an amplitude of approximately 80 milli-magnitudes. 
In addition to the eclipses of the close pair, the GIV compo- 
nent exhibits quasi-perio dic variability attributed to starspots 
Ivan Hamme et all 1 19941. 

A technique has been developed to obtain high preci- 
sion (10-20 //as) astrom etry of close stellar pairs (sep aration 
less than one arcsecond; Lane and Muterspaugh l2004ll ) using 
long-baseline infrared inter ferometry at t he Palomar Testbed 
Interferometer (PTI: IColavita et alJ I"l999ln . This level of pre- 
cision allows us to determine the orbital parameters of the 
V819 Her system at the level of ~ 3-5%, primarily lim- 
ited by the quality of available radial velocities. However, be- 
cause one can measure the astrometric motion of the close 
pair in the system, it is possible to directly determine the 
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relative inclination of the orbits. These measurements were 
made as part of the Palomar High-precision Astrometric Search 
for Exoplanet Systems (PHASES) program. PTI is located 
on Palomar Mountain near San Diego, CA. It was devel- 
oped by the Jet Propulsion Laboratory, California Institute of 
Technology for NASA, as a testbed for interferometric tech- 
niques applicable to the Keck Interferometer and other mis- 
sions such as the Space Interferometry Mission, SIM. It oper- 
ates in the J (1.2/mi), H (1.6^m), and K (2.2yum) bands, and 
combines starlight from two out of three available 40-cm aper- 
tures. The apertures form a triangle with one 110 and two 87 
meter baselines. 



2. Observations and Data Processing 

2.1. PHASES Observations 

V819 Herculis was observed using PTI on 31 nights 
in 2 003-200 5 using the observing mode described in 
lLane and Mutersp augh! l2004ll . This method for phase- 
referenced differential astrometry of subarcsecond binaries is 
briefly reviewed here. 

In an optical interferometer light is collected at two or more 
apertures and brought to a central location where the beams 
are combined and a fringe pattern produced. For a broadband 
source of central wavelength A the fringe pattern appears only 
when the optical paths through the arms of the interferometer 
are equalized to within a coherence length (A = A 2 /A/1). For 
a two-aperture interferometer, neglecting dispersion, the inten- 
sity measured at one of the combined beams is given by 



/ sin(7rx/A) 

I(x) = /o 1 + V — — sin (2nx/A) 

\ Tlx/ A. 



(1) 



where x is the differential amount of path between arms of the 
interferometer, V is the fringe contrast or "visibility", which 
can be related to the morphology of the source, and AA is the 
optical bandwidth of the interferometer assuming a flat optical 
bandpass (for PTI AA = 0.4yum). 

The location of the resulting interference fringes are related 
to the position of the target star and the observing geometry via 



d =lf ■ 5 s * +6 a (t,t) + c 



(2) 



where d is the optical path-length one must introduce between 
the two arms of the interferometer to find fringes. This quan- 
tity is often called the "delay." if is the baseline, the vector 
connecting the two apertures. is the unit vector in the source 
direction, and c is a constant additional scalar delay introduced 
by the instrument. The term 5 a fS \ ij is related to the differen- 
tial amount of path introduced by the atmosphere over each 
telescope due to variations in refractive index. For a 100-m 
baseline interferometer an astrometric precision of 10 /vas cor- 
responds to knowing d to 5 nm, a difficult but not impossible 
proposition for all terms except that related to the atmospheric 
delay. Atmospheric turbulence, which changes over distances 
of tens of centimeters, angles on order tens of arcseconds, and 
on subsecond timescales, forces one to use very short exposures 
(to maintain fringe contrast) and hence limits the sensitivity of 



the instrument. It also severely limits the astrometric accuracy 
of a simple interferometer, at least over large sky-angles. 

However, in narrow-angle astrometry one is concerned with 
a close pair of stars, and the observable is a differential astro- 
metric measurement, i.e. one is interested in knowing the angle 
between the two stars (A s = s| - ~s\). The atmospheric turbu- 
lence is correlated over small angles. If the measurements of 
the two stars are simultaneous, or nearly so, the atmospheric 
term subtracts out. Hence it is still possible to obtain high pre- 
cision "narrow-angle" astrometry IColavita|.ll994ll . 

To correct for time-dependent fluctuations in the atmo- 
spheric turbulence, observations consisted of operating PTI in 
a phase-referenced observing mode. After movable mirrors in 
the beam-combining lab apply delay compensation to the light 
collected from two 40 cm apertures, the light from each aper- 
ture is split using 30/70 beamsplitters. Seventy percent of the 
light is sent to the phase-tracking "primary" interferometric 
beam combiner which measures the time-dependent phase of 
one star's interferogram (fringes) caused by the atmospheric 
turbulence, and used in a feed-back loop to control the optical 
delay lines. 

The other 30% of the light is diverted to the "secondary" 
interferometric beam combiner. In this system there is an ad- 
ditional delay line with a travel of only * 500 microns. This 
is used to introduce delay with a sawtooth waveform with fre- 
quency on order a Hertz. This allows us to sample the interfer- 
ograms of all stars in the one arc second detector field whose 
projected separations are within the scan range. Laser metrol- 
ogy is used along all starlight paths between the 30/70 split 
and the point of interferometric combination to monitor inter- 
nal vibrations differential to the phase-referencing and scan- 
ning beam combiners. For V819 Herculis, the typical scanning 
rate in 2003 was one scan per second and four intensity mea- 
surements per ten milliseconds; these values were doubled in 
2004. The typical scan amplitude was 100 microns. An average 
of 3099 scans were collected each night the star was observed 
over a time span of 18 to 179 minutes. 

2.2. PHASES Data Reduction 

T he data reduction a lgorithm used was similar to that described 
in lLane and Muterspaugh! ll2004ll . but now all astrometric fit- 
ting is analyzed with differential right ascension and decli- 
nation as parameters, rather than projected separation. First, 
detector calibrations (gain, bias, and background) are applied 
to the intensity measurements. Next, a grid is constructed 
in differential right ascension and declination over which to 
search (in ICRS 2000.0 coordinates). For each point in the 
search grid we calculate the expected differential delay based 
on the interferometer location, baseline geometry, and time 
of observation for each scan. These conversions were sim- 
plified using the routines from the Naval Observatory Vector 
Astromet ry Subr o utines C Language Version 2.0 (NOVAS-C; 
see Kaplan et a l. 1 1989]). A model of a double-fringe packet 
is then calculated and compared to the observed scan to de- 
rive ax 2 value as a merit of goodness-of-fit; this is repeated for 
each scan, co-adding all of the x 2 values associated with that 
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point in the search grid. The model fringe template is found by 
observing single stars, incoherently averaging periodograms of 
their interferograms, and fitting a sum of Gaussians to the av- 
erage periodogram. This model effective bandpass is Fourier 
transformed into delay space to create a model interferogram. 
Sample data sets have been reanalyzed with a variety of model 
interferograms and the resulting astrometric solutions vary by 
less than one //as; this is largely due to the differential nature 
of the measurement. Note that in addition to the differential de- 
lay there are several additional parameters to the double fringe 
packet: fringe contrast and relative intensities as well as mean 
delay. These are all adjusted to minimize^ 2 on a scan-by-scan 
basis. The final x 2 surface as a function of differential right as- 
cension and declination is thus derived. The best-fit astrometric 
position is found at the minimum x 1 position, with uncertain- 
ties defined by the appropriate x 1 contour — which depends on 
the number of degrees of freedom in the problem and the value 
of the ^-minimum. 

One potential complication with fitting a fringe to the data 
is that there are many local minima spaced at multiples of 
the operating wavelength. If one were to fit a fringe model to 
each scan separately and average (or fit an astrometric model 
to) the resulting delays, one would be severely limited by this 
fringe ambiguity (for a 110-m baseline interferometer oper- 
ating at 2.2/mi, the resulting positional ambiguity is ~ 4.1 
milli-arcseconds). However, by using the ^-surface approach, 
and co-adding the probabilities associated with all possible de- 
lays for each scan, the ambiguity disappears. This is due to 
two things, the first being that co-adding simply improves the 
signal-to-noise ratio. Second, since the observations usually 
last for an hour or even longer, the associated baseline change 
due to Earth rotation also has the effect of "smearing" out all 
but the true global minimum. The final ^-surface does have 
dips separated by ~ 4.1 milli-arcseconds from the true loca- 
tion, but any data sets for which these show up at or above the 
4<x level are rejected. The final astrometry measurement and 
related uncertainties are derived by fitting only the Act region 
of the surface. 

The PHASES data reduction algorithm naturally accounts 
for contributions from photon and read-noise. Unmonitored 
phase noise shows up by increasing the minimum value of x 1 
surface. Comparison of this value with that expected from the 
number of degrees of freedom allows us to co-add the phase 
noise to the fit uncertainties. 

This method has been rigorously tested on both synthetic 
and real data. Data sets are divided into equal sized subsets 
which are analyzed separately. A Kolmogorov-Smirnov test 
shows the formal uncertainties from the PHASES data reduc- 
tion pipeline to be consistent with the scatter between subsets. 
After an astrometric solution has been determined, one can re- 
visit the individual scans and determine best-fit delay separa- 
tions on a scan-by-scan basis (the fringe ambiguity now be- 
ing removed). The differential delay residuals show normal 
(Gaussian) distribution, and Allan variances of delay residu- 
als agree with the performance levels of the formal uncertain- 
ties and show the residuals to be uncorrelated. It is concluded 
that the PHASES data reduction pipeline produces measure- 



ment uncertainties that are consistent with the internal scatter 
of the data on intra-night timescales. 

2.3. PHASES Astrometric Results 

The differential astrometry measurements are listed in Table 
^ in the ICRS 2000.0 reference frame. In order to evalu- 
ate the night-to-night astrometric repeatability of the data, the 
PHASES data were fit to a model consisting of a Keplerian 
orbit representing the Ba-Bb center of light (CL) motion and 
a low-order polynomial representing motion of the A-B orbit. 
The minimized value of reduced x 2 r = 4, implying either that 
the internal (i.e. derived from a single night of data) uncertainty 
estimates are too low by a factor of 2, or that the simple model 
is not appropriate for this system. Replacing the polynomial 
model for A-B with a Keplerian does not improve the value 
of xf- (to be expected given the limited fraction of the A-B or- 
bit covered by the PHASES data set). It is possible starspots 
caused astrometric jitter on this scale. The PHASES uncertain- 
ties presented in this paper have been increased by a factor 
of 2 to account for this discrepancy. The rescaled (raw) me- 
dian minor- and major-axis uncertainties are 15.2 (7.6) and 363 
(181) //as. The rescaled (raw) mean minor- and major-axis un- 
certainties are 19.6 (9.8) and 568 (284) //as respectively. 



2.4. Potential Systematic Errors 

The fractional precision of the PHASES astrometric measure- 
ments is ~ 10~ 4 ; at such an ambitious level there are many 
possible sources of systematic error that could appear on inter- 
night timescales. In particular, the system in question exhibits 
two potential astrophysical sources of measurement noise: 
starspots and eclipses. 



2.4.1. Starspots 

The m 40 milli-magnitude variability of V819 Herculis A has 
been attributed to starspots. Upper limits to shifts in the CL of 
a star caused by star spots are evaluated with a model com- 
prised of a uniform stellar disk (radius R) except for a zero- 
temperature (non-emitting) circular region of radius r tangent 
to the edge of the stellar disk (i.e. centered at x — R - r, y = 0). 
The CL is displaced by: 



c Lr I — *W* " k > f_ j—j^ ***** 



Rn(R 2 -r 2 ) 



r 2 /R 2 
1+r/R' 



(3) 



The non-emitting spots in this model would cause photo- 
metric variations proportional to the fractional area of the stel- 
lar disk covered: 



F r 
T ~ ~¥ 



(4) 



where F„ is the star's flux when no spots are present. Equations 
[3]and|3]provide a relationship between the apparent astrometric 
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TableQ 

PHASES data for V819 Herculis 



JD-2400000.5 


SRA 


SDec 




0"maj 


0e 


0RA 


°" Dec 


IT 1 

^"RA " Dec 


N 




(mas) 


(mas) 


(pas) 


(pas) 


(deg) 


(pas) 


(pas) 






53109.4777819 


49.6404 


-84.4971 


14.7 


566.0 


158.77 


527.6 


205.4 


-0.99704 


2011 


53110.4800679 


48.0940 


-84.1339 


23.7 


1199.8 


159.53 


1124.1 


420.3 


-0.99818 


1334 


53123.4555720 


49.1857 


-85.9322 


36.2 


1015.6 


162.47 


968.5 


307.9 


-0.99239 


1378 


53130.4397626 


48.4773 


-86.4138 


13.1 


411.8 


162.95 


393.7 


121.4 


-0.99359 


2537 


53137.4279917 


48.3926 


-87.1400 


27.9 


560.7 


164.34 


539.9 


153.7 


-0.98203 


1226 


53145.3928594 


48.3081 


-87.8018 


27.5 


633.0 


161.59 


600.7 


201.6 


-0.98962 


1673 


53168.3368802 


47.0273 


-89.7517 


30.1 


680.1 


162.93 


650.2 


201.7 


-0.98777 


1409 


53172.3496326 


47.3437 


-90.1341 


12.6 


339.9 


168.29 


332.9 


70.1 


-0.98310 


2560 


53173.3294522 


47.1580 


-90.3605 


16.0 


154.8 


33.97 


128.7 


87.5 


0.97549 


2904 


53181.3314386 


46.4330 


-90.7861 


15.0 


349.0 


169.71 


343.4 


64.1 


-0.97112 


2795 


53186.3020911 


45.6581 


-91.1217 


36.4 


853.0 


166.80 


830.5 


198.0 


-0.98202 


706 


53187.3022539 


46.1426 


-91.2229 


26.0 


882.1 


166.94 


859.3 


201.0 


-0.99114 


1578 


53197.2663645 


46.1848 


-92.1526 


9.5 


234.1 


164.87 


226.0 


61.8 


-0.98714 


5218 


53198.2404599 


46.2934 


-92.2558 


11.4 


109.3 


160.36 


103.0 


38.3 


-0.94842 


5404 


53199.2897673 


44.0250 


-91.9450 


49.2 


3283.0 


171.42 


3246.2 


492.3 


-0.99487 


946 


53208.2505295 


46.4281 


-92.4906 


13.2 


362.5 


37.67 


287.1 


221.8 


0.99719 


6558 


53214.2391404 


45.6333 


-93.3432 


10.9 


251.7 


169.45 


247.5 


47.3 


-0.97195 


5251 


53215.2293360 


45.6361 


-93.5176 


9.5 


221.3 


167.53 


216.1 


48.7 


-0.97963 


5723 


53221.2207072 


46.2536 


-92.9732 


17.6 


683.9 


38.91 


532.3 


429.8 


0.99861 


3998 


53228.2083438 


45.0879 


-94.3314 


14.4 


201.2 


169.45 


197.8 


39.5 


-0.92815 


3180 


53233.1820405 


45.0989 


-94.8462 


12.0 


129.5 


167.67 


126.5 


30.0 


-0.91190 


3303 


53234.2006462 


44.8264 


-94.7640 


15.2 


75.7 


172.74 


75.1 


17.9 


-0.51361 


3701 


53235.2168202 


45.2148 


-94.9186 


17.0 


214.3 


176.57 


213.9 


21.2 


-0.60018 


2094 


53236.1665478 


44.4594 


-94.8865 


9.5 


156.1 


166.59 


151.9 


37.4 


-0.96553 


6684 


53481.5043302 


30.4429 


-103.3356 


22.2 


622.2 


38.18 


489.3 


385.0 


0.99730 


3301 



Table 1. PHASES data for V819 Herculis. All quantities are in the ICRS 2000.0 reference frame. The uncertainty values 
presented in this data have all been scaled by a factor of 2 over the formal (internal) uncertainties within each given night. 
Column 6, <f> e , is the angle between the major axis of the uncertainty ellipse and the right ascension axis, measured from increasing 
differential right ascension through increasing differential declination (the position angle of the uncertainty ellipse's orientation 
is 90 - <p e ). The last column is the number of scans taken during a given night. The quadrant was chosen such that the larger 
fringe contrast is designated the primary (contrast is a combination of source luminosity and interferometric visibility). 



and photometric shifts caused by star spots. For stars of typical 
radius 1 milli-arcsecond, the simplified model gives a roughly 
linear relationship of 0.8 pas of maximum astrometric shift per 
milli-magnitude of photometric variability. 

The largest possible astrometric shift by a star spot is given 
by evaluating a slightly different model. In this case, the star 
spot fills the (non-circular) area from the star's edge to a chord 
at distance x from the star's true center. The astrometric shift 
is 



~R 



3(l + arcsinf + f (l - f ) ) 
with corresponding photometric variations of 



(5) 



F 



n . x 

— h arcsin — 
2 R 



x 
~R 



R 2 



2\l/2\ 



(6) 



In the V819 Herculis system, the effect of a single, cold 
starspot is maximally ~ 25pas assuming a stellar radius of 0.8 
mas. If multiple starspots cause the variability, this effect is re- 
duced. In addition, the effect of limb-darkening is to further 
reduce the astrometric error. 



2.4.2. PHASES Observations During Ba-Bb Eclipses 

Using the published size s and temperatures for Ba and Bb from 
Ivan Hamme et alJ Il994ll . it is found that the magnitude of the 
astrometric shift in CL position during eclipse compared to 
what it would be outside of eclipse can be greater than 100 
pas. Because this shift is larger than PHASES astrometric mea- 
surement precisions, six measurements taken during eclipse are 
omitted from the data tables and the fits. 



2.5. Previous Differential Astrometry Measurements 

Previously published differential astrometry measur ements of 
the A- B (wide) system have been tabulated by Hartkopf et alJ 
[2004] in the Fourth Catalog of Interferometric Measurements 
of Binary Stars. In several cases discrepancies were found be- 
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Fig. 1. The maximum effect of star spots on astrometric mea- 
surements versus the photometric variations they cause. 



tween uncertainties quoted in the original works (or uncertainty 
estimates omitted in the catalog); in these cases the uncertainty 
estimates from the original works are used. All of these mea- 
surements were made using the technique of speckle interfer- 
ometry. These measurements are included in the combined fit 
to help complete coverage of the A-B orbit. 

Many of the previous differential astrometry measurements 
were published without any associated uncertainties. To allow 
these to be used in combined fits with other data sets, the aver- 
age uncertainties were determined as follows. The uncertainties 
were initially assigned values of 10 milli-arcseconds in sep- 
aration and 1 degree in position angle. A (single) Keplerian 
model was fit to the data, and residuals in separation and po- 
sition angle were treated individually to remove outliers and 
update the uncertainty estimates. This procedure was iterated 
until uncertainties were found consistent with the scatter. A 
double Keplerian model (as in eq. [8] to allow for the Ba-Bb 
subsystem) does not improve the fit; the measurements are in- 
sensitive to this small signal. These 22 data points have average 
uncertainties of 5.92 milli-arcseconds and 0.689 degrees. 

A Keplerian model was fit to the data points for which un- 
certainty estimates were available to determine whether these 
were systematically too large or too small, and to find outliers. 
The uncertainty estimates are found to be systematically too 
small; this factor was larger in position angle than in separa- 
tion. Upon iteration, it was found that the separation uncertain- 
ties for these 12 data points needed to be increased by a factor 
of 2.43 and the position angle uncertainties by 3.53. Again no 
improvement was seen in fitting to a double Keplerian model. 

These previously published measurements are listed in 
Tables|2]and|3] 



2.6. Radial Velocity Data 

A large number of radial velocity measurements of components 
A and Ba from four observatories were reported in lScarfe et alJ 
ll994ll . lScarfeetal] indicate several measurements as outliers; 
these measurements have not been used. There are 72 com- 
ponent A velocities and 50 component Ba velocities in the 
McDonald/Kitt Peak data set, 70 component A and 49 com- 
ponent Ba velocities in the DAO set, and 92 component A and 
90 component Ba velocities in the DDO data set. 

The velocity measurements for each of three data sets pre- 
sented were fit to double Keplerian models separately to de- 
termine the average velocity uncertainties (measurements from 
McDonald and Kitt Peak were mixed together in the original 
w ork, and these were analyzed together as one group). As noted 
in lScarfe et al.llll994ll . the velocity precisions for component A 
differed from those of Ba. Uncertainties were derived for each 
data set by fitting to a double Keplerian model and examining 
the scatter in the residuals for each component separately. The 
uncertainty guesses were updated and the procedure iterated. 
The initial values for the component A velocities were 0.46 
kms~' for the McDonald/Kitt Peak and DAO data sets, and 
0.92 km s for the DDO set; for component Ba, all were set to 
2 kms -1 . The average A and Ba uncertainties are 0.43 kms~' 
and 1.955 km s _1 for the McDonald/Kitt Peak velocities, 0.465 
km s _1 and 3.025 km s for the DAO measurements, and 1.015 
km s _1 and 3. 105 km s _I for the DDO observations. 



3. Orbital Models 

Basic models have been applied to the astrometric data. The 
simplifying assumption was made that the Ba-Bb subsystem 
is unperturbed by star A over the timescale of the observing 
program, allowing the model to be split into a wide (slow) in- 
teraction between star A and the center of mass (CM) of B, and 
the close (fast) interaction between stars Ba and Bb. The results 
presented in this paper result from modeling both the A-B and 
Ba-Bb motions with Keplerian orbits. 

In general, one cannot simply superimpose the results of 
the two orbits. The observable in the PHASES measurements 
is the separation of star A and the CL of the Ba-Bb subsystem. 
Because the CL of Ba-Bb, the CM of Ba-Bb, and the loca- 
tion of star Ba are generally all unequal, a coupling amplitude 
must be added to the combined model. This coupling ampli- 
tude measures the relative size of the semi-major axis of the 
Ba-Bb subsystem to that of the motion of the CL of the Ba-Bb 
subsystem. The sign of the superposition is determined by the 
relative sizes of the mass and luminosity ratios of the stars Ba 
and Bb. As an example, if the CL is located between the CM 
of Ba-Bb and the location of star Ba, the motion of the CL will 
be in opposite direction to the vector pointing from Ba to Bb. 
For a subsystem with mass ratio M B b/M Ba and luminosity ratio 
^Bb/^Ba, the observed quantity is 



Jobs = 7"a-B - 



M Bb /M Ba - L Bb /L Ba 

(1 + M Bb /M Ba )(l + L Bb /L Ba ) 



'"Ba-Bb 



(7) 



where r^-s is the model separation pointing from star A to the 
CM of B, and r Ba -Bb is the model separation pointing from star 
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Tabled 

Previous astrometry with uncertainties 



Year 


p 


9 






Reference 




(mas) 


(deg) 


(mas) 


(deg) 




1980.4766 


36 


278.1 


4.86 


7.06 


McAlister et al. [1983] 


1980.4820 


47 


265.7 


4.86 


7.06 


McAlister et al. H9831 


1982.2924 


107 


154.7 


2.43 


1.76 


Fuetal. T19971 


1991.3192 


46 


252.4 


7.29 


0.71 


Hartkopf et al. r 19941 


1992.3104 


48 


123.4 


7.29 


0.71 


Hartkopf et al. r 19941 


1992.6086 


69 


141.2 


7.29 


7.06 


Balega et al. 119991 


1993.2058 


96 


149.9 


7.29 


0.71 


Hartkopf et al. r 19941 


1994.7078 


102 


169.8 


1.49 


1.06 


Hartkopf et al. [2000] 


1995.3138 


86 


177.9 


7.29 


0.71 


Hartkopf et al. ri9971 


1999.5681 


115 


164.0 


6.90 


3.21 


Horch et al. r20021 


1999.5681 


109 


160.4 


6.90 


3.21 


Horch et al. r20021 


2000.5563 


101 


175.6 


6.90 


3.21 


Horch et al. [2002] 



Table 2. Previous differential astrometry data with published uncertainties for V8 19 Herculis. p uncertainties have been increased 
by a factor of 2.43 and those for 9 by a factor of 3.53. In many cases 9 has been changed by 180 degrees from the value appearing 
in the original works. 



Tabled 

Previous astrometry without uncertainties 



Year p 9 Reference 
(mas) (deg) 


Year p 9 Reference 
(mas) (deg) 


1981.4568 0.065 134.1 McAlister et al. [1984] 


1984.7009 0.076 188.3 McAlister et al. [1987a] 


1981.4678 0.061 134.2 McAlister et al. ri9841 


1985.4816 0.051 224.4 McAlister et al. r 1987a! 


1981.4706 0.063 135.2 McAlister et al. ri9841 


1985.5228 0.047 228.3 McAlister et al. r 1987b! 


1981.4732 0.062 135.6 McAlister et al. ri9841 


1986.4079 0.034 354.9 Balega et al. ri9891 


1982.5027 0.109 155.2 McAlister et al. 11987al 


1986.6454 0.038 109.0 Baleeaet al. H9891 


1983.0702 0.112 163.0 McAlister et al. 11987al 


1987.2673 0.080 140.9 McAlister et al. [19891 


1983.4202 0.105 165.8 McAlister et al. H987al 


1988.2529 0.107 157.5 McAlister et al. fl989! 


1983.7151 0.107 171.3 McAlister et al. H987al 


1988.6599 0.108 162.4 McAlister et al. r 19901 


1984.3732 0.100 180.2 McAlister et al. H987al 


1989.2305 0.101 170.1 McAlister et al. r 19901 


1984.3760 0.088 181.4 McAlister et al. 11987a! 


1989.7058 0.091 177.7 Hartkopf et al. 11992! 


1984.3840 0.087 181.5 McAlister et al. 11987a! 


1990.2651 0.074 190.2 Hartkopf et al. [19921 



Table 3. Previous interferometric differential astrometry data without published uncertainties for V819 Herculis. The uncertain- 
ties presented in this table were determined by the scatter in the data. For the values in this table, all uncertainties were taken to 
be cr p = 5.92 mas and erg = 0.689 degrees. 



Ba to star Bb. Including this coupling term for astrometric data 
is important when a full analysis including radial velocity data 
is made. 

Alternatively, one can directly combine a model of the A-B 
system with a model of the motion of the CL of Ba-Bb. For 
purely astrometric data such a model is appropriate. In this 
case, there is no sign change for the Ba-Bb CL model, and no 
extra coupling amplitude is required. This solution can be used 
to provide a high-precision apparent orbit which can be used 
to compare different instruments. This model differs from the 
previous one in that the true semi-major axis of the Ba-Bb orbit 
is not a parameter, but is replaced by the apparent photocentric 
motion directly. 



fobs = rA-B + '"Ba-Bb.C.O.L. (8) 



4. Orbital Solution and Derived Quantities 



The best-fit combined astrometry-radial velocity orbital solu- 
tion produces a set of parameters listed in Table@] The reduced 
X 2 r of the combined fit to PHASES, radial velocity, and previ- 
ous differential astrometry data is 1.33. This combined set has 
521 degrees of freedom with 20 parameters. This value for^^ 
is slightly higher than one would expect, but this is likely due to 
the manner in which the uncertainties had to be derived. All pa- 
rameter uncertainties have been increased by a factor of Vl.33 
to reflect this difference. It is noted that the value for the in- 
clination of the Ba-Bb pair of 79.0 ± 3.3 degrees agrees with 
that determined by eclipse observations (8 1 .00 + 0.36 for the 
sy nchronous and 80.63 + .33 for the asynchronous solutions 
of Ivan Hamme et al.l Il994ll ). Also presented is a table of de- 
rived parameters of direct astrophysical interest (0. A fit to the 
astrometric data alone does not constrain several of the orbital 
parameters, but does constrain apparent semi-major axes of the 
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TableEl 

Orbital models for V819 Herculis 



Table© 



PHASES 
+ Pre. + RV 


Pa-b (days) 


2019.79 ±0.36 


T 0A -b (MJD) 


52628.1 ±1.3 


eA-B 


0.6731 ±0.0015 


Ia-b (degrees) 


57.09 ±0.22 


u>a-b (degrees) 


42.55 ±0.23 


Q A -b (degrees) 


322.40 ±0.16 


P Ba -Bb (days) 


2.2296337 ±1.9 X lO" 6 


T ,Ba-Bb (MJD) 


52627.18 ±0.30 


e Ba-Bb 


0.0041 ±0.0033 


iBa-Bb (degrees) 


79.0 ±3.3 


0)Ba-Bb (degrees) 


47 ±48 


&Ba~Bb (degrees) 


312.9 ±4.8 


V ,m/k (kms^ 1 ) 


-3.388 ±0.059 


V ,dao (km s" 1 ) 


-3.373 ±0.064 


V ,ddo (km s~') 


-3.35 ±0.12 


M A (M ) 


1.765 ±0.095 


M Ba+ Bb (M G ) 


2.512 ±0.067 


M B b/M Ba 


0.757 ±0.020 


l^Bbl ' Lsa 


0.261 ±0.045 


d (parsecs) 


67.96 ±0.87 



Table 4. Orbit models for V819 Herculis. Pre.: Previous differ- 
ential astrometry measurements, listed in Tables |5] and [3] The 
luminosity ratio LsblLBa is for K-band observations. 

Table© 

Derived physical parameters 



Parameter 


Derived Value 


a A -B (AU) 


5.076 ± 0.046 


"Ba-Bb (AU) 


0.04541 ± 0.00040 


cos O (degrees) 


23.6 ±4.9 


M Ba (M Q ) 


1.430 ±0.041 


M Bb (M Q ) 


1.082 ±0.033 



Table 5. Physical parameters for V819 Herculis derived from 
the combined orbital solution. 



wide and narrow orbits to good precision: oa-b = 73.9 + 0.6 
mas and asa-Bb = 108 + 8/ias (Ba-Bb CL orbit). 

The combined fit A-Ba-Bb orbit is plotted is Figure |5J 
PHASES measurements of the Ba-Bb orbit CL motion is plot- 
ted in Figure [3] Residuals to the combined fit are shown in 
Figures |4] [5] and |6j no evidence for additional system com- 
ponents is observed. 

4.1. Mutual Inclination 

The mutual inclination O of two orbits is given by 

cos <£ = cos z'i cos 12 + sin ii sin cos (Q; - Q.2) (9) 

where i\ and z'2 are the orbital inclinations and fii and Q.2 are 
the longitudes of the ascending nodes. From the combined or- 
bit solution a value for this angle of 23.6 + 4.9 degrees is de- 
rived for the V819 Her system. This low value is below the 
limit for inclination-eccentricity oscillations derived by Kozai 



Known Mutual Inclinations 


Star 


Mutual Inclination 


Reference 




(degrees) 




V819Her 


23.6 ±4.9 


This paper 


k Peg 


43.8 ± 3.0 


[Mutersoauah et al., 2005] 


r] Vir 


30.8 ± 1.3 


rHummel et al.. 20031 


6 Hya 


39.4 


THeintz. 19961 


£UMa 


132.1 


THeintz. 19961 


Algol 


98.8 ±4.9 


TLestrade et al.. 19931 






[Pan et al., 1993] 



Table 6. Unambiguously known mutual inclinations of triple 
systems. The value for Algol i s determined usin g the measure- 
ment precisions and values of P an et alJ lll993ll for all but the 
A-B nodal position angle of 52 ± 5 degrees from lLestrade et all 
ll 993ft 



(1962; 39.2 degrees), and is consistent with the small measured 
value of the eccentricity of the Ba-Bb pair. 

The mutual inclination of the orbits of triple systems is 
of particular interest for und erstanding the conditions un - 
der which the system formed ISterzik and Tokoviniri Eo021. 
Without both radial velocity and visual (or astrometric) orbits 
for both systems in a triple, unambiguous determinations of the 
longitudes of the ascending nodes (and thus of the mutual in- 
clination) are impossible. To date there has only been a very 
small number of cases where the mutual inclination can be un- 
ambiguously determined (Tabled 

With the tally of systems for which unambiguous mutual 
inclinations have been determined now at six, it is reasonable 
to consider the distribution of t hese orientations. The previou s 
work on this subject is that of Sterzik and Tokovinin 12002], 
who determined theoretical distributions resulting from a va- 
riety of initial conditions within star forming regions. At the 
time, the authors cited only three systems from which mutual 
inclin ations were know n, listed in a previous work by one of 
them I Tokoviniri Jl993|. (One o f these three, £ Cnc, is most 
recently listed by Heintz 1 1996] as still having an ambiguous 
mutual in clination, and is not incl uded here.) For comparison to 
real stars. ISterzik and Tokovinin! instead included the 22 t riple 
systems in the Multiple Star Catalog iTokovininl [l997ll for 
which both visual orbits were known, but the ascending nodes 
had 180 degree ambiguities. To correct for this lack of informa- 
tion, for each system they included both possible mutual incli- 
nations in a combined cumulative distribution (this distribution 
is referred to as ST). This procedure assumes the ambiguity is 
divided evenly between the lower and higher possible angles; 
i.e. that an equal number of the "true" mutual inclinations are 
the lesser of the two possible angles as are the greater. 

In Figure the continuous distribution function of mutual 
inclinations for the six unambiguous systems (this distribution 
is referred to as 6U) is plotted with ST and the theoretical distri- 
bution for random orientations (referred to simply as Random). 

The two-sided Kolmogorov-Smirnov (K-S) probability for 
agreement between 6U and ST is 0.46; the one-sided KS 
probability between ST and Random is 0.07 and that for 6U 
and Random is 0.04. Thus, the 6U set confirms the result of 
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Previous Astrometry PHASES Astrometry 




Orbit Phase (Degrees) Orbit Phase (Degrees) 



Fig. 2. The orbit of V819 Herculis. (Top left) Previous (speckle) differential astrometry measurements with derived uncertainty 
ellipses. (Top right) One season of PHASES astrometry. (Bottom left) CM velocities of the wide pair. (Bottom right) Radial 
velocity of the Ba component. The A-B motion has been removed for clarity. 



ISterzik and Tokovininl i2002ll that mutual inclinations are not 
consistent with random orientations and show a slight prefer- 
ence for coplanarity. The sets 6U and ST agree much better 
with each other than either do with random orientations, but 
the agreement probability is still low; this is likely due to the 
assumption inherent to forming the ST set by including both 
possible orientations, which dilutes the distribution away from 
coplanarity. A greater number of systems is required to better 
determine the distributions, and observational selection effects 
should also be considered. 

5. Conclusions 

PHASES interferometric astrometry has been used to measure 
the orbital parameters of the triple star system V819 Herculis, 
and in particular to resolve the apparent orbital motion of the 
close Ba-Bb pair. The amplitude of the Ba-Bb CL motion is 



only 108 ± 9 //as, indicating the level of astrometric precision 
attainable with interferometric astrometry. By measuring both 
orbits one is able to determine the mutual inclination of the 
two orbits, which is found to be 23.6 ± 4.9 degrees. Such a low 
mutual inclination implies a lack of Kozai oscillations. 

Further improvement in determining the system distance 
and component masses will require improved radial velocity 
data. Given that the A component is evolved and the Ba-Bb sys- 
tem undergoes eclipses and hence the B components have ac- 
curately measured radii, this system may become a very fruitful 
laboratory for high-precision testing of stellar models. 
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axis (left) and along an axis rotated 165 degrees East of North (equivalent to position angle 285 degrees; right); this is the median 
position angle of the minor axis of the PHASES uncertainty ellipses. Because the orientation of the PHASES uncertainty ellipses 
varies from night to night, no single axis is ideal for exhibiting the PHASES precisions, but this median axis is best aligned to do 
so. The wide A-B orbital motion has been removed for both plots. The error bars plotted have been stretched by a factor of 2 over 
the formal uncertainties as discussed in the text. The high ellipticity of the uncertainty ellipses causes neither the right ascension 
nor the declination uncertainties to be near the precision of the minor axis uncertainties, which have median uncertainty of 15.2 
//as. For clarity, measurements with projected uncertainties larger than 200 //as are not shown in the plots. Zero phase indicates 
periastron passage. 
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Fig. 4. Residuals for PHASES differential astrometry of V819 Herculis. The error bars plotted have been stretched by a factor 
of 2 over the formal uncertainties as discussed in the text. The high ellipticity of the uncertainty ellipses causes neither the right 
ascension nor the declination uncertainties to be near the precision of the minor axis uncertainties, which have median uncertainty 
of 15.2 /j.as. Due to the roughly North-South alignment of the baseline used for most of the measurements, the more sensitive axis 
was typically declination. The right ascension and declination plots show only those points for which the projected error bar is less 
than 1 milli-arcsecond. The bottom left plot shows the residuals along a direction that is 165 degrees from increasing differential 
right ascension through increasing differential declination (equivalent to position angle 285 degrees), which corresponds to the 
median direction of the minor axis of the PHASES uncertainty ellipses; only measurements with uncertainties less than 500 /ias 
along this axis are plotted. Because the orientation of the PHASES uncertainty ellipses varies from night to night, no single axis 
is ideal for exhibiting the PHASES precisions, but this median axis is best aligned to do so. The bottom right plot shows residuals 
along the minor axis of each measurement's uncertainty ellipse. 
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Fig. 5. Residuals to the combined model for previous differential astrometry of V8 19 Herculis. 
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Fig. 7. Cumulative distribution of the observed distribu- 
tion of angles between angular momentum vectors of the 
six systems for which unambiguous mutual inclinations have 
been determi ned. This is c ompared with the results from 
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distribution. Also shown is the theoretical distribution for ran- 
dom orientations. 
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